Pourquoi cette modélisation ? Pourquoi la Rochelle ?

Au départ : articuler la distribution spatiale des résidents et des opportunités

• notamment les emplois, mais aussi les écoles, les commerces, les espaces verts

• comprendre les prix de l’immobilier et la différence qu’induit l’accessibilité à ces opportunités

Pour aller plus avant, les émissions de CO₂ par exemple, il faut passer du potentiel à des flux réalisés (en espérance). C’est cette question qui s’est posée lors d’un travail pour l’agglomération de la Rochelle (sur le périmètre du SCoT La Rochelle Aunis).

3 documents de travail

Parodi et Timbeau (2023) : description du modèle théorique, discussion du modèle gravitaire, analyse des propriétés de MEAPS

Parodi et Timbeau (2024b) : estimation de MEAPS et comparaison avec le modèle gravitaire. discussion de la stratégie d’estimation

Parodi et Timbeau (2024a) : analyse du lien entre densité et émission et construction de modèles pour les fréquences, les détours et les modes de transport

un site https://preview.meaps.fr sur lequel vous trouverez les documents de travail et les liens vers les dépôts .

Accessibilité aux emplois

A la Rochelle: accessibilité aux emplois par commune

Ce qui demande de localiser emplois et résidents

• Individus : données carroyées INSEE c200m 2017 et 2019

• Emplois : MOBPRO (recensement fichier détail) 2019 imputés par commune et par secteur NAF 5 en fonction des surfaces professionnelles issues des fichiers fonciers et du RFP

De l’accessibilité aux émissions de CO₂ (km parcourus)

Comment produire cette carte ?

Principe général de la modélisation

Du tenseur \([km^m_{ijk}]\) on peut alors calculer l’espérance des kilomètres par carreau \(km_{ij}\), pour une catégorie particulière \(km_{k}\), pour un mode spécifique \(km^k\), etc…

Synthèse : Pourquoi MEAPS est mieux

• Théoriquement MEAPS décrit un processus d’appariement spatialisé explicite et vraisemblable. La saturation permet un appariement complet [ce qui n’est pas le cas pour gravitaire].

• Empiriquement, il faut utiliser un modèle de Poisson ou une fonction objectif entropie relative de Kullback-Leibler pour estimer un modèle de flux (MEAPS ou gravitaire).

• Empiriquement, MEAPS 0p. peut servir de modèle de référence. En introduisant des odds ratios, on peut paramètrer MEAPS pour représenter une préférence locale (par ex.).

• Empiriquement, le modèle gravitaire quant à lui repose sur des effets fixes ou une procédure de Furness pour fonctionner. Mais on triche avec les masses alors que les masses sont l’argument principal.

• Empiriquement, avec seulement 2 ou 3 paramètres, on arrive à expliquer (presque) aussi bien avec MEAPS qu’avec un gravitaire à effets fixes ou Furness

• Hors échantillon, on peut reproduire les flux avec seulement la localisation de l’emploi et des résidents

• En ajoutant une information infra-communale, on améliore la capacité de MEAPS à repduire les flux, mieux qu’un gravitaire saturé.

MEAPS : principes

On part d’une analogie radiative plutôt que gravitaire : le milieu traversé compte (stouffer1940? ; simini2012? ; Simini, Maritan, et Néda 2013)

1. absorption : Chaque individu part de \(i\) et rencontre les opportunités \(j\) classées dans l’ordre des distances. A chaque rencontre il a une probabilité \(p_i\)d’être absorbé. Le voisinage compte, puisque si je traverse un milieu dense en opportunités je vais moins loin

2. saturation : Chaque fois qu’un individu est absorbé, l’opportunité est diminuée. Il y a concurrence entre les individus pour accéder aux opportunités. A nouveau le voisinage compte : si je suis dans un milieu dense en individus je dois aller plus loin

3. priorité : la saturation est résolue en donnant la priorité au premier arrivé. On donne un ordre aux individus et on détermine leur absorption en fonction de 1. et de 2. dans l’ordre de priorité

4. ergodicité : on tire un grand nombre d’ordres aléatoirement, en répétant 1 2 et 3 pour obtenir en moyenne une allocation indépendante de l’ordre arbitraire, ce qui achève la prise en compte du voisinage

Par construction ce modèle respecte le principe de séparabilité. La saturation permet en outre d’assurer que chaque individu occupe un emploi et que chaque emploi est occupé si les opportunités sont toutes accessibles à tous les individus.

Radiatif plutôt que gravitaire : des arguments théoriques

Le modèle gravitaire dans sa forme générale :

\[ f_{ij} =c\times \frac{n_i^\alpha \times e_j ^\beta }{d_{ij}^\delta } \]

est un modèle qui ne tient pas compte du voisinage : le flux est généré par une masse au départ, une autre à l’arrivée mais l’entourage n’intervient pas (fotheringham1983? ; simini2012? ; stouffer1940?).

On peut définir une notion de séparabilité qui n’est respectée à l’origine que si \(\alpha=1\) et à la destination si \(\beta=1\) : si on divise un groupe au même point \(ij\) en deux sous groupes, on veut que la somme de flux générés par les deux sous groupes soit exactement le flux généré par le groupe entier si ils ont les mêmes paramètres (de comportement). De la même façon, si on sépare un groupe d’individus ou d’emplois très proches en des sous groupes à une résolution plus élevée on veut que le changement de résolution n’est qu’un impact faible sur les flux modélisés (à la limite nul).

Radiatif ou gravitaire : Poisson versus log normal

On estime le modèle gravitaire par cette équation en supposant un processus générateur log normal Lenormand, Bassolas, et Ramasco (2016) :

\[ log(f_{ij}) = \alpha \times log(n_i) + \beta \times log(e_j) - \delta \times d_{ij} + c + \varepsilon_{ij} \;\;;\; \;\varepsilon_{ij} \sim \mathcal{N}(0, \sigma^2) \]

C’est une mauvaise façon de faire (Flowerdew et Aitkin 1982) parce que le processus générateur est celui d’un Poisson (comptage) :

\[ P(\hat f_{ij}= f_{ij}) = \frac{e^{-\hat f_{ij}}\times \hat f_{ij}^{f_{ij}}}{f_{ij}!} \;et\; log(\hat f_{ij}) = \alpha \times log(n_i) + \beta \times log(e_j) - \delta \times d_{ij} +c \]

ou, la comparaison de deux densités doit se faire par l’entropie relative de Kullback-Leibler (kullback1951?). La log vraisemblance d’une table de contingence s’écrit :

\[ \mathcal L = \sum_{ij} f_{ij} log(\hat f_{ij} /n) \] au lieu de l’ajustement usuel par les MCO (on peut s’en approcher par des MCO pondérés par \(f_{ij}\) (Agresti 2002, pp.146‑148).

Radiatif ou gravitaire : effets fixes ou aléatoires

Reste que les estimations ne s’arrêtent pas là : on dispose généralement d’une information que l’on veut utiliser (ou respecter). On connait les \(n_i\) et le \(e_i\). L’estimation par Poisson ou par résidu log normaux ne permet pas le respect de cette information. On ajoute donc généralement des effets fixes ou aléatoires (suivant l’hypothèse que l’on pose pour la projection par exemple) :

\[log(\hat f_{ij}) = \alpha \times log(n_i) + \beta \times log(e_j) - \delta \times log(d_{ij}) + log(a_i) + log(b_j) \]

Dans ce cas le modèle, ne permet pas d’identifier \(\alpha\) ou \(\beta\) – on les pose à 1 pour respecter le principe de séparabilité.

On peut également utiliser une procédure de Furness (itérative, puis estimation non linéaire, on montre que \(\alpha=\beta=1\) ) :

\[ a_i = \frac{n_i}{\sum_j {{b_j n_i ^ \alpha e_j ^ \beta} / { d_{ij}^{\delta} }}} = \frac{n_i ^ {1-\alpha}}{\sum_j { b_j e_j ^ \beta / d_{ij}^{\delta} }} \]

\[ b_j = \frac{e_j}{\sum_i { {a_i n_i ^ \alpha e_j ^ \beta} / { d_{ij}^{\delta} }}} = \frac{e_j ^ {1-\beta}}{\sum_i {a_i n_i ^ \alpha / d_{ij}^{\delta} }} \]

Radiatif ou gravitaire : ajustement de MEAPS

MEAPS est un modèle sans paramètre (comme simini2012?). On peut ajouter des paramètres dans le modèle :

\[ \tilde{p}_{abs,ij} = \frac{c_{abs} \times \omicron_{ij}} {1+c_{abs} \times \omicron_{ij}} \]

Puis en définissant une structure \(O\) sur les \(\omicron_{ij}\) si on ne veut pas d’un modèle saturé (\(\mathcal{L}\) est bien sûr l’entropie relative KL ;)).

\[ \hat \theta = \underset{\theta}{\mathrm{argmin}} \, \mathcal{L}(f^{meaps}_{ij}(O(d_{ij},\theta))) \]

Estimations et ajustement sur MOBPRO : une surprise !

En intégrant les informations infracommunales, en simulant MEAPS, en agrégeant au niveau communal, on produit une meilleure estimation :

Détours

Détours : Boucles simples

\[ logit(Prob^{simple}_k) = \alpha + \beta \times typmen_k +\varepsilon_k \\ avec \varepsilon_k \sim Binomial \]

En utilisant L’EMP 2019, on repère les boucles simples (trajet domicile travail puis retour).

Détour : Boucles complexes

Pour une distance \(d\) domicile travail, la longueur de la boucle est \(L = K \times d\). En posant \(\gamma\) la proportion de détours rapportée à la distance domicile travail on a \(\gamma = K - 2\). On fait une régression par quantile :

\[ log(\gamma_l) = \alpha + \beta \times log(1 + d_l) + \varepsilon_l \\ avec\ \varepsilon_l \sim \mathcal{N} \\ en\ minimisant \sum_l \left\lvert{\varepsilon_l}\right\lvert \]

En utilisant L’EMP 2019, on repère pour les boucles complexes où le premier ou le dernier trajet est entre le domicile et le travail la longueur de la boucle et la distance domicile travail (boucles B ou D). On mesure alors \(\gamma\).

Nombre de boucles

On modélise les boucles par un modèle de Poisson augmenté pour la sous-dispersion Conway-Maxwell-Poisson (Conway1962?)

\[ P(x=n|\lambda,\nu) = \frac{\lambda^n}{Z(\lambda, \nu) \times (n!)^\nu} \quad pour\ n=0,1,2... \]

avec l’équation suivante :

\[ log(nb^{bcl}_k) = \alpha + \gamma \times log(L^{bcl}_k) + \delta \times dens_{res, k} + \tau \times voiture_k + \varepsilon_k \\avec \\ \varepsilon_k \sim CMP(\lambda, \nu) \]

Fréquences : Nombre de boucles

En utlisant l’EMP 2019, on compte par individu le nombre de boucles.

Modes de transport

On suit (mcfadden1974d?) en estimant un Random Utility Model

\[ U_l^m = \alpha_m + \beta \times tt_{bcl,m} + \gamma_m \times log(L_{bcl}) + \delta_m \times log(NdV_l) \\ \lambda_m \times dens_l + \mu_m \times typmen_l + \nu_m \times voiture_l + \varepsilon_{m,l} \\ avec\ \varepsilon_{m,l} \sim Loi\ de\ Gumbel \]

La source est l’EMP 2019. On n’utilïse pas le recensement qui interroge sur le mode habituel [on pourrait vérifier quand même].

Carte au carreau 200m

Lien entre niveau de vie et distance parcourue

Scénarios de densification : le lien entre la densité et la forme urbaine

La densité et les émissions de CO2

En regardant finement la distribution spatiale des individus on s’affranchit du problème de l’autosélection.

En résumé

On combine des sources de données :

• La localisation des résidents au carreau 200m (données carroyées INSEE 2019)

• La localisation des emplois au carreau 200m (recensement fichier détail MOBPRO + fichiers fonciers + RFP, 5 secteurs)

• Les réseaux de transport pour calculer les distances et les temps de parcours par mode entre chaque paire

• les flux inter communaux (recensement fichier détail MOBPRO) pour ajuster le modèle de flux infra communal

• l’enquête EMP 2019 pour évaluer la longueur des boucles, leur complexité et leur fréquence

• l’enquête EMP 2019 pour évaluer les modes

On modélise les flux, les fréquences, les détours, les modes sans agréger les individus spatialement, avec des propriétés pertinentes et explicites.

Sur ces sources ouvertes (ou presque), nous proposons une approche en open source afin d’assurer la transparence et la reproductibilité.

On obtient un modèle calibré qui peut être utilisé avec une certaine confiance pour intrapoler, interpoler, analyser la géographie (localisations, réseaux) et faire des analyses de scénarios

Extensions possibles

• Utiliser les données enrichies par le SDES sur les émissions dans l’EMP 2019

• Introduire des dépendances entre les termes de l’équation fondamentale : par exemple, les flux différenciés par ménages, par catégories d’emploi, par mode

• Inclure la congestion dans les calculs de temps de parcours en voiture (par ex. données Mapbox sur le trafic usuel, données de flotte potentiellement) ou le temps de transport effectif moyen (GTFS temps réel disponible pour certaines agglomérations).

• D’autres structures dans les odds de MEAPS, des processus stochastiques non binomiaux

• Affiner les catégories de ménages, celles d’emploi

• Exploiter la tension sur l’emploi

Futurs développements: Marseille et le passage à l’échelle

Actuellement, nous travaillons sur la métropole d’Aix-Marseille-Provence. C’est un passage à l’échelle significatif : La matrice de flux de la Rochelle est 5000 \(\times\) 5000 \(\approx\) 25M. Celle de Marseille est de 26000 \(\times\) 38000 \(\approx\) 1G, soit 40 fois plus grande. On calcule par exemple les paires origines destination à une vitesse entre 5000 et 1000 paires/s/vCPU.

Futurs développements: Autres motifs

Un des enjeux est calculer l’espérance pour les autres motifs que le motif professionnel.

• mobilités scolaires : recensement fichier détail plus ensemble d’opportunité plus simple, pas de saturation

• mobilités pour le commerce : pas de source origine destination, ensemble d’opportunités très large, très divers. Nous explorons l’EMC² AMP pour une estimation en construisant une couche commerce synthétique (un motif agrégé, pondéré par les usages) et en ayant des matrices O/D partielles.

• mobilités pour les autres motifs (santé, démarches, loisir, socialisation) : pas de source origine destination, si ce n’est L’EMC². Motifs encore plus divers que pour les commerces. Exploration d’une couche synthétique construite sur la population avec une estimation des matrices O/D partielles.

Futurs développements: Autres données

• L’utilisation d’autres données comme celles de (lévy2023?) (bornage de téléphone, qui permettent des matrices o/d à la résolution de l’IRIS). Un projet ANR est en cours (de candidature) pour collaborer avec Lévy et Coldefy sur cette question (Digitalisation et décarbonation des mobilités - MOBIDEC).

• La future chaire Dynamiques Urbaines a pour objet d’accueillir ces travaux et de centraliser les ressources (dont les données de mobilité).

• D’autres traces numériques peut éventuellement permettre d’affiner un ou plusieurs des modèles employés (notamment MEAPS) et de les confronter à des résolutions infra-communales.

• L’utilisation de Fidéli sur le CASD permet d’avoir une information fine (revenu, structure familiale au carreau 200m) qui fournirait une source importante pour la projection et la pondération des kilomètres. Il y a également la possibilité d’un début d’analyse temporelle. Mais les règles actuelles de préservation de l’anonymat empêchent de produire et de publier une information dérivée au carreau 200m à partir de ces données.